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Abstract 

We consider a simple model of lossless interaction between a two-level single atom and 
a standing-wave single-mode laser field which creates a one-dimensional optical lattice. 
Internal dynamics of the atom is governed by the laser field which is treated to be classical 
with a large number of photons. Center-of-mass classical atomic motion is governed by 
the optical potential and the internal atomic degree of freedom. The resulting Hamilton- 
Schrodinger equations of motion are a five-dimensional nonlinear dynamical system with 
two integrals of motion, the total atomic energy and the length of the Bloch vector are 
conserved during the interaction. In our previous papers the motion of the atom has 
been shown to be regular or chaotic (in the sense of exponential sensitivity to small 
variations in initial conditions and/or the system's control parameters) in dependence 
on values of the control parameters, the atom-field detuning and recoil frequency. At 
exact atom-field resonance, exact solutions for both the external and internal atomic 
degrees of freedom can be derived. The center-of-mass motion does not depend in this 
case on the internal variables, whereas the Rabi oscillations of the atomic inversion is a 
frequency modulated signal with the frequency to be defined by the atomic position in the 
optical lattice. We study analytically correlations between the Rabi oscillations and the 
center-of-mass motion in two limiting cases of a regular motion out off the resonance: (1) 
far-detuned atoms and (2) fastly moving atoms. The main focus of the paper is chaotic 
atomic motion that may be quantified strictly by positive values of the maximal Lyapunov 
exponent. It is shown that atom, depending on the value of its total energy, can either 
oscillate chaotically in a well of the optical potential or fly ballistically with weak chaotic 
oscillations of its momentum or wander in the optical lattice changing the direction of 
motion in a chaotic way. In the regime of chaotic wandering atomic motion is shown to 
have fractal properties. We find a useful tool to visualize complicated atomic motion - 
Poincare mapping of atomic trajectories in an effective three-dimensional phase space 
onto planes of atomic internal variables and momentum. The Poincare mappings are 
constructing using a translational invariance of the standing laser wave. We find common 
features with typical non-hyperbolic Hamiltonian systems — chains of resonant islands of 
different sizes imbedded in a stochastic sea, stochastic layers, bifurcations, and so on. The 
phenomenon of sticking of atomic trajectories to boundaries of regular islands, that should 
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have a great influence to atomic transport in optical lattices, is found and demonstrated 
numerically. 

PACS 42.50.Vk, 05.45.Mt, 05.45.Xt 

1 Introduction 

Light exerts mechanical forces on matter. This hypothesis was suggested by Kepler pQ in 1619 
to explain a deviation of the comet's tails flying nearby the Sun. It was Maxwell who in 1873 
estimated the light pressure, using his theory of electromagnetism [2], and has shown that 
it is very small. Peter Lebedev was the first who in 1899 measured the light pressure on a 
macroscopic body [3]. The first experiments on deviation of microscopic particles by light have 
been carried out by W. Gerlach and O. Stern [3], by P. Kapitza and P. Dirac and by O. 
Frisch jHj. 

Manipulation of atomic motion with the help of laser beams, creating an optical lattice, 
is one of the most fastly growing field of modern physics (for a review see, for example, [Zj). 
There are different theoretical and experimental aspects of this interaction including cooling and 
trapping of atoms, Bose-Einstein condensation, quantum computing and processing information 
with atoms. 

In this paper we review our recent results on nonlinear coherent dynamics of a single two- 
level atom in an optical lattice created in a one- dimensional cavity by two counterpropagating 
laser waves. We are working in the strong-coupling regime and neglect all the losses. We show 
that even in a one-dimensional approximation the atomic motion can be very complicated. 
We analyze both regular and chaotic motion of atoms in a stationary standing-wave laser field 
containing a large number of photons. It should be stressed that there is a difference between 
various types of erratic atomic motion in an optical lattice. Chaotic motion is strictly defined as 
a motion of a deterministic nonlinear dynamical system that is exponentially sensitive to small 
variations in the system's initial conditions or/and its control parameters. There are different 
types of chaotic motion of atoms in an optical lattice, including chaotic nonlinear oscillations 
of atomic center of mass in a well of the optical potential, chaotic ballistic motion, when the 
atomic momentum oscillates chaotically around a value of the average momentum, and the last 
but not least, chaotic wandering of an atom when it changes its direction of motion in a chaotic 
way jHJ H3 El EEU EH EH]- All the types of chaotic atomic motion are quantified by positive 
values of the maximal Lyapunov exponent. 

In an optical lattice chaotic motion in the strict sense of this notation may occur when 
there is no any kind of noise, including atomic spontaneous emission which is a random process. 
The respective deterministic atomic equations of motion are approximate ones, but they are 
fundamental since spontaneous emission may be considered as a quantum noise. Random 
walking is a kind of motion that occurs with ultracold atoms which are detuned far away from 
the carrier laser frequency, so their internal degrees of freedom can be eliminated adiabatically. 
Because the values of the momentum of ultracold atoms are compared with the value of the 
photon momentum, each time after emitting a spontaneous photon atom gets a kick in a random 
direction. This effect is a quantum analogue of the classical random walking. 

In a general situation we should take into account both the internal atomic motion and 
spontaneous emission events. In this case, however, the equations of motion cease to be a 
deterministic dynamical system because they include random terms and one may expect much 
more complicated type of atomic motion which, besides of chaotic motion, caused by the fun- 
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damental atom-field interaction, includes a purely stochastic component caused by random 
events of spontaneous emission. We have shown recently that in a range of the control param- 
eters (detuning, laser intensity, and recoil frequency) and initial conditions atoms may change 
their direction of motion erratically even if their momenta are much larger than the photon 
momentum. We will call this type of motion as chaotic walking. 

The main aim of this paper is to describe different aspects of deterministic atomic motion 
in an optical lattice, both regular and chaotic ones. The effects of spontaneous emission on the 
atomic motion will be considered in a forthcoming paper. 



2 Hamilton- Schrodinger equations of motion 

We consider a two-level atom with mass m a and transition frequency u a , moving with the 
momentum P along the axis X in an ideal cavity through the standing laser wave with the 
field frequency ujf and the wave vector kf. In the frame, rotating with the frequency ujf, the 
standard cavity QED Hamiltonian is the following: 

P 2 1 

H = 1 — h(uj a — ujf)a z — Ml (er_ + <r + ) cos kfX. (1) 

2777" d 2 

Here a± >z are the Pauli operators which describe the transitions between lower, |1), and upper, 
1 2), states. fl is the Rabi frequency which is proportional to the square root of the number of 
photons in the wave y/n. The standing- wave field must be strong enough (n ^> 1), so we can 
neglect a back reaction of atoms on it and consider the field classically. For electronic degree 
of freedom the simple wavefunction is 

m)) = <w)+KW), (2) 

where a and b are the complex- valued probability amplitudes to find the atom in the states |2) 
and |1), respectively. Using the Hamiltonian we get the Schrodinger equation 



.da uj a -Uf 

i— = a — izocos kfX, 

Ojv Zi 

.db ujf-u a 

i— = — o — I la cos ktX, 

dt 2 1 ' 



(3) 



where the atomic position X is considered as a parameter. Let us introduce instead of the 
complex-valued probability amplitudes a and b the following new real-valued variables: 

u = 2Re(ab*), v = -21m(ab*) , z = \a\ 2 - \b\ 2 , (4) 

which are the quadratures of the atomic dipole moment (u and v) and the atomic population 
inversion, z. 

In the process of emitting and absorbing photons, atoms not only change their internal 
electronic states but their external translational states change as well due to the photon recoil. 
If the atomic average momentum is large as compared to the photon momentum hkf, one 
can describe the translational degree of freedom classically satisfying to classical Hamilton 
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equations of motion. The dynamics in the strong-coupling regime is now governed by the 
Hamilton-Schrodinger equations 



x = u r p, 

p = —u sinx, 

u = Av, (5) 
v = —Au + 2z cos x, 
z = — 2-ucosx, 

where x = kf(X) andp = (P)/hkf are classical atomic center-of-mass position and momentum, 
respectively. Dot denotes differentiation with respect to dimensionless time r = Qt. The 
normalized recoil frequency, uj t = hk 2 /m a Q <C 1, and the atom-field detuning, A = (uf—u a )/fl, 
are the control parameters. The system has two integrals of motion, namely the total energy 

W = ^p 2 + U, (6) 

where 

A 

U = — ucosa: — —z, (7) 

is the potential energy, and the Bloch vector 

u 2 + v 2 + z 2 = l. (8) 

The conservation of the Bloch vector length immediately follows from Eqs. (@J). 

Equations © with two integrals of motion constitute a Hamiltonian autonomous system 
with two degrees of freedom and motion on a three-dimensional hypersurface with a given energy 
value W . Generally, such a system has a positive Lyapunov exponent A, a negative exponent 
equal in magnitude to positive one, and two zero exponents. The sum of all Lyapunov exponents 
of a Hamiltonian system is zero . The maximal Lyapunov exponent characterizes the mean 
rate of the exponential divergence of initially close trajectories, 

A=limA(r), A(r) = lim -In ^4, (9) 

r->oo v ; ' K ' <5(0H0T 5(0)' V ' 

and serves as a quantitative measure of dynamical chaos in the system. Here, 5{r) is a distance 
(in the Euclidian sense) at time r between two trajectories close to each other at initial time 
t = 0. The dependence of A on control parameters has been calculated in [HUH] with the similar 
system. It has been shown that dynamical chaos in a strongly-coupled atom-field system exist 
in a wide range of parameters and initial atomic momentum p . The result of computation of 
the maximal Lyapunov exponent with our system © in the space of control parameters, uo r 
and A, is shown in Fig.0 In white regions of the map maximal Lyapunov exponent A is almost 
zero and the dynamics is stable. In other regions the positive values of A show the Lyapunov 
instability. 

In all the numerical simulations we use physically realistic parameters and initial conditions. 
For example, we can chose cesium atoms with the transition wavelength A a ~ 852 nm. The 
Rabi frequency Q depends on many factors including the field strength, which could be changed 
in a wide range. In most computations we shall use the Rabi frequency of Q ~ 10 GHz. With 
this value of the Rabi frequency we get the normalized recoil frequency to be uj r = 10~ 5 . Also 
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Figure 1: Maximal Lyapunov exponent A vs the detuning A and the normalized recoil frequency 
u r : po = 200, z = — 1, Uq = Vq = 0. 

we put the initial position xq = 0. The detuning A could be varied in a wide range, and the 
Bloch variables are restricted by the length of the Bloch vector (jHJ). The most interesting effects 
are observed with rather cold atoms. For example, po = 200 taken in computing A in Fig. ^ 
with our normalization corresponds to the atomic velocity v a pa 0.7 m/s. It should be noted 
that we use in this paper the normalization to the laser Rabi frequency Q, not to the vacuum 
(or single-photon) Rabi frequency as it has been done in our previous papers [HI HH E31 EB EG] • 
So the ranges of the normalized control parameters, taken in this paper, differ from those in 
the cited papers. 



3 Regular dynamical regime 
3.1 Exact atom- field resonance 

At exact resonance, A = 0, one can easily find an additional integral of motion, 

u = const = uq. (10) 

In this case the fast and slow variables are separated from each other allowing one to integrate 
exactly the reduced equations of motion. Total energy becomes equal to 

Wr = yp 2 - M COSX, (11) 

and the potential energy gets the simple form 

Ur = —Mo cos x. (12) 
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The center-of-mass translational motion of the atom in such a spatially periodic potential of 
the standing wave is described by the simple nonlinear equation for a free physical pendulum 



x + u r u sin x — 0, 



(13) 



and does not depend on evolution of the internal degrees of freedom. 

The translational motion is trivial when uq is zero. In spite of the zero potential field, a 
structure of a standing wave is still present in a cavity. In this case, the atom will move in one 
direction with a constant velocity, and the Rabi oscillations modulated by the standing wave 
will occur. In general case one can easily get from (|TT| the dependence p(x) 



V 



— (Wr + uo cosx) 



(14) 



which gives the phase portrait of the system in the plane (x,p) (Fig. EK)- It is the phase 
portrait of a nonlinear pendulum with three types of trajectories depending on the value of its 
energy Wr: oscillator-like motion in a potential well if Wr < uq, a separatrix if Wr = uq, and 
ballistic-like motion if Wr > uq. 



(a) 



(b) 



P 





Figure 2: The typical regular phase portraits for the translational degree of freedom: (a) 
A = 0; (b) |A| >0.2. 



For the initial values xq = and xq = u> r po, the equation for the translational motion (fT3*|) 
has the solution 

2 arcsin [K sn [^cu r UoT, K]] , K 2 ^ 1; 



x(t) 



2 am 



Po dn 



(15) 



K 2 > 1, 



Po cn [y/u r u T, K], K 2 ^ 1; 



(16) 



K 2 >1, 



where 



2 V u o 



is the modulus of the elliptic Jacobi functions. The solution gives the critical value of the 
atomic momentum 

Per = 2^U /U r . (18) 

Atoms with po ^ Per are trapped by the standing-wave field, the result that is well-known 
from early studies The modulus K is simply connected with the normalized value of the 
difference between the energy of the atom and its value on the separatrix 

K* = 1+ Wr - Uo . (19) 
2u 

As to internal atomic evolution, it depends on the translational degree of freedom since the 
force of the atom-field coupling depends on the position of the atom in a periodic standing-wave 
potential. The equation for the atomic population inversion z(r) is derived from the two last 
equations of the set © with A = 0: 



z = T-2yi - z 2 -u§cos[z(t)], (20) 

where cos[x(t)] is a known function of the translational variables only which can be found with 
the help of the exact solutions obtained. It is easy to find the exact solution of Eq. (J2H|) 

z ( r ) — T\j l-^o sin ^2 J cos xdr' + tp^j , (21) 
where the sign is opposite to that for the initial value Zq and 

ipo = =F arcsin Z ° (22) 
V 1 ~ u o 

is an integration constant. The internal energy of the atom could be considered as a frequency- 
modulated signal with the instant frequency 2cos[x(r)] and the modulation frequency x = 
u} r p(r), but it is correct only if the first value is much greater than the second, i. e. \w r po\ ^ 2. 
Such a signal is shown in Fig. for a ballistic atom (p = 5000, v a ps 17.5 m/s). With 
|avpo| ^ 2 the modulation disappears and the signal becomes a periodic one with the frequency 
ijj r p. With fast atoms, |covp | ^ 2 and p ~ p ^> p cr (Raman-Nath approximation) Eq. (j2~Tj) is 



simplified 



z(t) ~ z smuo r p T ^ sm oo r Po T - (2oJ 

VrPo ^rPo 



3.2 Non-resonant Rabi oscillations 

With comparatively small detunings A the dynamics of slow atoms can be chaotic. In this case 
the Rabi oscillations are still a signal with a frequency modulation, but the amplitude is not 
constant anymore, it jumps chaotically with the characteristic time l/u r p. With comparatively 
large detunings, the Rabi oscillations become regular, but with a prominent periodic amplitude 
modulation, while the frequency modulation is not so deep (Fig. Eb). 
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Figure 3: Rabi oscillations: (a) frequency modulation at exact resonance, A = 0, po = 5000; 
(b) amplitude modulation far from resonance, A = —4, po = 5000; (c) Doppler-Rabi resonance, 
A = -4, po = 400000. In all panels, z = u = v / 0~l>, v = 0. 



In two limit cases, |A| ^> max[|u; r p|, 2] and \oJ r p\ 3> max[|A|,2], the analytic solutions can 
be obtained. In this both cases u and v are the harmonic functions with the frequency A, and 
for the internal energy we have 



z + 



z 



2u 2^/u 2 + 



A 



cosx sin(Ar + <po), |A| 3> max[|a; r p|, 2], 



2y^ 



(24) 



o + ^o 



u r p 



cos(Ar + <p ) sinu; r p r, |^ r p | ^ max[|A|, 2], uJ T p\ 4, 



where ip = arcsin(w 



u : 



+ Vq). The Rabi oscillations now are amplitude-modulated signals 



with two characteristic frequencies \uJ r p\ and |A|. The larger frequency is the main frequency 
and the other one is the modulation frequency. In the solution for fast atoms we also used the 
Raman-Nath approximation, x ~ u r poT, that is correct if the initial kinetic energy uo T p^/2 is 
much greater than the amplitude of potential energy variations (equal to 2 in our case). With 
A = 0, the solution has the form (}2*3*j) . but without the last term which is small. Solutions 
(}2"4*|) show a good correspondence with the numerical experiments, performed in ^2] with the 
similar equations. The typical amplitude-modulated Rabi oscillations are shown in Fig. Eb- 

More exact solution for u can be found using the approximation z ~ const ~ Zq (i. e. 
z max — z min |^o|); which is correct in the both limit cases considered above, excluding z ~ 0. 
Then from Eq. (jSJ) we get the equation for a driven linear oscillator 



u + A 2 u ~ 2z A cos x, 



(25) 



S 



which has the solution 

u(t) ~ 2z sin At J cos At cos xc/r — 2z cos Ar y sin Ar cos x<ir+ 

+«o cos Ar + f o sin At. (26) 

For |A| 3> \w r p\, the solution can be approximated as follows: 

2zq I 

u w cosx + yu 2 , + Vq sin(Ar + <p ). (27) 

Using Eqs. (j21| and (|2*7|). we get the periodic potential with the spatial period 7r (in difference 
from the resonant potential with the period 2ir): 

2zq 

U ~ — — - cos 2 x + const. (28) 

The corresponding phase portrait is shown in Fig. |2b- 

When the frequencies are close, \u r p\ ~ |A|, the Doppler-Rabi resonance takes place ^2] i n 
spite of the fact that the detuning may be very large. Let us consider the standing wave as a 
combination of two counter-propagating waves. In the frame, moving with the atomic velocity, 
their frequencies, uj\ and u 2 , are different because of the Doppler effect: 

UJ\ = LUf LUf, UJ2 = OJf H OJf, (29) 

where v a is the atomic velocity and c is the speed of light. The atom is rather slow so we 
can neglect the relativistic effects. Let us consider atoms fast enough for the Raman-Nath 
approximation p pa p to be valid. Renormalizing all the frequencies to Q, we define the 
dimensionless detunings between the atomic transition and the running wave frequencies as: 

Ai = — - — = A-uj r p , A 2 = — - — = A + LO r p . (30) 

The condition |A| = |w r po| leads to the resonance between the atom and one of the waves. If 
|A| ^> 1, we can neglect the interaction with the other wave and consider the atom as if only 
one wave with the frequency uj\ or 0J2 exists. In the field of the wave, say, cui, the dynamics can 
be described by the Bloch-like equations 

ii = Aiv, v = — Aiu + z, z = —v, (31) 

in which the interaction energy does not depend on the atomic position and its amplitude value 
is twice smaller as compared to the standing wave. Eqs. (j31|) have the solution 

m Ai v ( A? 1 \ 
z = — 5-(l - coster) sinc^r + z — H n cosu z t , (32) 

where u z = a/A 2 + 1 = a/(A — uj r po) 2 + 1. At the exact Doppler resonance (Ai = 0), atomic 
internal energy z oscillates with the dimensionless frequency 1, and the amplitude of oscillations 
is maximal. Numerical simulations with Eqs. (0) shows that this speculations are correct 
(Fig. |2b, where po = 400000 and v a ~ 1400 m/s), and even very far from the resonance A = 
the deep Rabi oscillations can be observed for the atoms to be fast enough. 
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Figure 4: A typical atomic trajectory in the regime of chaotic wandering: xq = 0, po = 300, 
zq = — 1, uq = v = 0, u r = 10~ 5 , A = —0.05. 

4 Irregular dynamics: chaos and fractals 
4.1 Chaotic atomic wandering 

In Fig. ^ we depict the maximal Lyapunov exponent map in the space of control parameters, 
uj r and A. The maximal Lyapunov exponent depends not only on the parameters uo r and A, 
but on initial conditions of the system (jSJ), as well. Especially important is a value of the initial 
momentum, p . The most interesting effects can be observed with rather cold atoms, when the 
initial atomic kinetic energy is close to the amplitude of the optical potential. In this case we 
get the chaotic wandering of an atom in the standing wave. A typical chaotic atomic trajectory 
is shown in Fig. EJ 

It follows from (jSJ) that the translational motion of the atom at A ^ is described by the 
equation of a nonlinear physical pendulum with a frequency modulation 

x + uj r u{r) sinx = 0, (33) 

in which u is the function of all the other dynamical variables. The normalized Rabi oscillation 
frequency is a value of the order of uj' z = y/ A 2 + 4 which substantially exceeds the frequency 
of small- amplitude translational motion y/u r uo <C 1 in the potential well. Taking this into 
consideration, the mechanism of the arising of chaos can be revealed ^Oj- The stochastic layer 
width was estimated as 

— © 3 -6£>- (34) 

where uj = ^2oo r \A\/uj' z , uj' z /ujq ^> 1. The D value is the energy change in the neighbourhood 
of the unperturbed separatrix normalized with respect to the pendulum separatrix energy uo^. 
Small changes in the energy causes comparatively small changes in the frequency of oscillations. 
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For the energies of motion that are strongly different from the separatrix energy, that is, close to 
potential well bottoms and high above optical potential U hills, small frequency changes cause 
small phase changes during the translation motion period. However, close to the unperturbed 
separatrix, where the period of oscillations tends to infinity, even small frequency changes can 
cause substantial phase changes. This is the reason for the exponential instability of motion 
of the parametric nonlinear oscillator (|33J) and chaotic atomic motion in the field of a periodic 
standing wave. 

A clear idea of the character of chaotic wandering can be developed using the model of "two 
potentials". At resonance, the optical potential U reproduces the structure of the standing 
wave in the cavity (jl2j) with the 2ir period (the phase portrait in Fig. Far from the 

resonance, the potential has the period n and is approximately described by Eq. (J25J) and the 
corresponding phase portrait is shown in Fig. Eb- These potentials will be called resonant and 
nonresonant, respectively. We can say that, when the motion in the cavity is chaotic, the both 
potentials "virtually" coexist. The well depths in both structures change as time passes, and 
an atom randomly gets into one or another structure every time when it crosses a standing 
wave node. The probability of getting into the resonant or nonresonant potentials depends on 
the detuning. Near the resonance atom is in the resonant potential almost the whole time and 
only rarely gets into nonresonant one for a short time. 

In our study ^2] we have shown that chaotic wandering has fractal properties. 

4.2 Dynamical atomic fractals 

In Fig. El we depict the scheme of a gedanken experiment that consists of a Fabry-Perot optical 
microcavity with two detectors and cold atoms to be placed in the cavity. To avoid complica- 



atom 

standing-wave field 




X 



mirrors 



detectors 



Figure 5: The schematic diagram shows a standing- wave microcavity with detectors. 

tions that are not essential to the main theme of this work, we consider the cavity with only two 
standing-wave lengths. Atoms, one by one, are placed at the point xq = with different values 
of the detuning A. We measure a time when an atom reaches one of the detectors, the exit 
time T, and study the dependence T(A) under the other equal conditions imposed on the atom 
and the cavity field. Fig. |U] shows the function T(A) for atoms with the initial momentum 
Po = 200 (v a pa 0.7 m/s). The exit time function demonstrates an intermittency of smooth 
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Figure 6: Atomic fractals. Exit time of cold atoms T vs the detuning A: po = 200, zq = — 1, 

Uq = V = 0. 
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curves and complicated structures that cannot be resolved in principle, no matter how large 
the magnification factor. The middle and low panels in Fig. El show successive magnifications 
of the function for the small intervals. Further magnifications reveals a self-similar fractal-like 
structure that is typical for Hamiltonian systems with chaotic scattering [To] . 

The exit time T, corresponding to both smooth and unresolved A intervals, increases with 
increasing the magnification factor. Theoretically, there exist atoms never reaching the detec- 
tors inspite of the fact that they have no obvious energy restrictions to leave the cavity. Tiny 
interplay between chaotic external and internal dynamics prevents these atoms from leaving 
the cavity. The similar phenomenon in chaotic scattering is known as dynamical trapping. In 
[T2"] for the similar fractal we have computed the Hausdorff dimension and shown that it is not 
integer. 

Different kinds of atomic trajectories before detection can be characterized by the number 
m — 1 of changing the sign of momentum. An m-th trajectory corresponds to the atom which 
changes the direction of motion before being detected m — 1 times. There are also special 
separatrix-like trajectories following which atoms asymptotically approach to the points with 
the maximum of potential energy, having no more kinetic energy to overcome it. In difference 
from the separatrix motion in the resonant system (A = 0) with the initial atomic momentum 
p cr , a detuned atom can asymptotically reach one of the stationary points even after several 
oscillations in a well. Let us define the m.S-trajectory as a trajectory when the atom changes the 
direction of motion m — 1 times and then begin the separatrix-like motion. Such asymptotical 
motion takes the infinite time, so the atom will never be detected. 

The smooth A intervals in the first-order structure (Fig.® upper panel) correspond to atoms 
which never changes the direction of motion, i. e. m = 1, and reaching the right detector. The 
unresolved singular points in the first-order structure with T = oo at the border between 
the smooth and unresolved A intervals are generated by the lS'-trajectories. Analogously, 
the smooth and unresolved A intervals in the second-order structure (Fig. |UJ middle panel) 
correspond to the 2-nd order and the other trajectories, respectively, with singular points 
between them corresponding to the 2S'-trajectories and so on. 

There are two different mechanisms of generation of infinite detection times, namely, dy- 
namical trapping with infinite oscillations (m = oo) in a cavity and the separatrix-like motion 
(m ^ oo). The set of all detunings generating the separatrix-like trajectories is a countable 
fractal. Each point in the set can be specified as a vector in a Hilbert space with m integer 
nonzero components. One is able to prescribe to any unresolved interval of m-th order struc- 
ture a set with m integers, where the first integer is a number of a second-order structure to 
which trajectory under consideration belongs in the first-order structure, the second integer 
is a number of a third-order structure in the second-order structure mentioned above, and so 
on. Such a number set is analogous to a directory tree address: " <a subdirectory of the root 
directory>/<a subdirectory of the 2-nd level>/<a subdirectory of the 3-rd level>/...". Un- 
like the separatrix fractal, the set of all detunungs leading to dynamically trapped atoms with 
m = oo seems to be uncountable. 

The scattering function in the regime of chaotic wandering, time of exit T, depends in a 
complicated way not only on the control parameters but initial conditions as well. In Fig. [7|we 
demonstarte the view of this function, whose values are modulated by color, in two coordinates, 
the initial atomic momentum p and the atom-fireld detuning A. From the fragment (a) to 
the fragment (f ) we increase subsequently the resolution. One can see increasing complexity of 
the scattreing function with a prominent self-similarity. The computation has been performed 
with the recoil frequency uj r = 9.17 • 10~ 5 . 
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4.3 Poincare mapping 

The five variables in the equations of motion (jHJ) minus the two integrals of motion (JOJ) and 
(JHJ) provide motion in a three-dimensional space. To visualize the motion we use the idea of 
mapping trajectories on two-dimensional planes. Since we have no time-periodic perturbations 
in our equations of motion (0) we cannot map trajectories through equal intervals of time 
provided by a period of a perturbation. However, the system has a characteristic space period 
27r imposed by the standing wave. So we map trajectories on a chosen plane at those time 
moments when atoms reach the positions when cosx = 1. We close our phase space along the 
position variable with the period 2ir. The condition cos x = 1 under fixed values of the integrals 
of motion © and (jHJ) defines a closed two-dimensional surface in the phase space points of which 
characterize unambigiously the system's states. In other words, there is a set of points on this 
two-dimensional surface which corresponds to each trajectory with a given value of the energy 
W. This set can be projected onto a plane of any system's variables except for the position 
x. Such a projection is, generally speaking, two-valued because the two-dimensional surface is 
closed. However, one can map trajectories in its "eastern" and "western" parts separately. 

In Fig. |Hk and b we demonstrate the Poincare mappings of a number of atomic trajectories 
in the western (u < 0) and eastern (u > 0) hemispheres of the Bloch sphere (u, v,z) on the 
plane v — z, respectively. We fix the values of the detuning A = —0.05, the recoil frequency 
uj r = 10 -5 , the total energy W = 33.8, the initial position x = 0, and map the trajectories 
with different other initial conditions, which are restricted by © and (JHJ). All the mappings 
were obtained with ballistic atoms whose momenta slightly (but chaotically for some initial 
conditions) oscillate around the effective value p c s = 2600 that corresponds to the chosen value 
of the energy. It is such a value of the momentum which an atom has at the moments when 
its potential energy U is zero. In Fig. El we demonstrate the Poincare mappings of a number of 
atomic trajectories in the western (u < 0) and eastern (u > 0) hemispheres of the Bloch sphere 
(u, v, z) on the plane v — z just like as in Fig. [HI with W = 33.8 but with another value of the 
total energy W = 36.45 and the effective momentum p e g = 2700. A series of bifurcation occurs 
just between these values of energy and we get in the end a central critical point instead of a 
saddle. One can see a typical structure with surviving nonlinear resonances of different orders 
around the center point and overlapping resonances. 

In Figs. |H] and El a general views of the mappings in the western and eastern hemispheres are 
shown. The pictures are rather typical with chaotic Hamiltonian systems [TJj . We see regions of 
regular motion in the form of islands and chains of islands filled by regular trajectories which are 
known as Kolmogorov-Arnold-Moser (KAM) invariant curves. The islands are imbedded into 
a stochastic sea, and they are produced by nonlinear resonances of different orders. Increasing 
the resolution of the mapping, one can see that big islands are surrounded by islands of a 
smaller size each of which, in turn, is surrounded by a chain of even more smaller islands, and 
so on to infinity. Stochastic layers of the oo-like form are situated between the islands. From 
the physical point of view, they are formed by broken and overlapping nonlinear resonances. 
From the mathematical point of view, a stochastic layer is a heteroclinic structure formed by 
transversal intersections of stable and unstable manifolds of hyperbolic stationary points. A 
fractal-like structure of generations of islands, a trademark of Hamiltonian chaos, is clearly 
seen on projections of motion in both the western and eastern hemispheres. To illustrate what 
happens under increasing the resolution of the Poincare mapping, we plot in Fig. |St a zoom of 
a small region in the stochastic layer in Fig. |Hk. 

We would like to pay attention to another typical phenomena in Hamiltonian systems, so 
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W = 33.8 



u > 



z o 




0.65 - 



Figure 8: Poincare mapping in the Bloch variable space, (a) u < (western Bloch hemisphere), 
(b) u > (eastern Bloch hemisphere), (c) magnification of the small region in (a) fragment, 
(d) mapping with a single chaotic trajectory in (b) fragment, illustrating the effect of sticking: 
W = 33.8, p eff = 2600, u r = 10" 5 , A = -0.05. 
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u < W = 36.45 u > 




v v 



Figure 9: Poincare mapping in the Bloch variable space, (a) u < (western Bloch hemisphere), 
(b) u > (eastern Bloch hemisphere). The parameters are the same as in Fig- El but W = 36.45, 
Pes = 2700. 

called sticking fTJ EH E3 EH1 EI] I n Fig- Etl we demonstrate the phenomenon of sticking 
in the eastern Bloch hemisphere. The trajectory shown demonstrates an intermittent type of 
motion. It wanders for a while in the stochastic chaotic trajectory, whose instability 

is characterized by a positive value of the finite-time maximal Lyapunov exponent. Then it is 
sticked to the boundaries of the outmost visible chain of regular islands, where it is practically 
regular with zero value of the respective finite-time maximal Lyapunov exponent. It may take 
a large amount of time to find a gap in a cantori structure, surrounding the outmost KAM 
tori, and to get off in the stochastic sea. The process is repeated as time grows. It should be 
stressed that sticking influences strongly transport properties in Hamiltonian systems giving 
rise to anomalous diffusion, algebraic tails in distributions of the Poincare recurrence times and 
of times and lengths of atomic flights. 

In Fig. ITUk and b we map the same atomic trajectories as in Figs. |H1 and H3 onto the plane 
p—z. In this case both the parts of the closed two-dimensional surface have the same projections 
because the set is symmetric with respect to the hyperplane v = 0. 

In order to quantify instability of the trajectories on the Poincare mappings in Figs. IHland 
we have computed the maps of the maximal Lyapunov exponents exactly with the same 
initial conditions and parameters as in those figures. The results in the vq — zq coordinates for 
eastern hemispheres at W — 33.8 and W = 36.45 with u > are shown in Figs. ITTk and b, 
respectively. A rather good correspondence between the Poincare mapping and the maximal 
Lyapunov exponents proves that the Poincare mapping we have constructed is a good means to 
visualize compicated dynamics of the coupled internal and external atomic degrees of freedom. 

The Poincare mapping with a rich structure of regular and chaotic regions is typical only 
in a range of values of the total energy W. At the values W > 40 atoms move regularly and 
the respective Poincare mapping consists of regular invariant curves only. With decreasing 
the energy, a series of bifurcations occurs with appearing of resonant islands of different order. 
When decreasing the energy even more, global chaoticity takes place. At exact resonance A = 0, 
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W = 33.8 



W = 36.45 





Figure 11: The maximal Lyapunov exponent vs the initial values of the Bloch variables Vq and 
z with u > 0. (a) W = 33.8, p cS = 2600, (b) W = 36.45, p oS = 2700. The other values are 
the same as in Fig. |H1 
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the value W = Uq corresponds to a separatrix in mechanical variables. Out off resonance this 
separatrix is broken, and atoms may wander chaotically in the optical lattice with a respective 
irregular Poincare mapping. At W < u + Az /2 (including negative values) atoms are trapped 
in wells of the optical potential and oscillate there. 

5 Conclusion 

We have considered a simple model of lossless interaction between a two-level single atom and a 
standing-wave single-mode laser field which creates a one-dimensional optical lattice. Analytical 
solutions of the Hamilton-Schrodinger equations of motion have been derived and analyzed in 
some limiting cases of regular atomic motion. Correlations between the Rabi oscillations and 
the center-of-mass motion have been established and demonstrated. In the regime of chaotic 
wandering the atomic motion has been shown to have fractal properties. Using a special type 
of the Poincare mapping of atomic trajectories in an effective three-dimensional phase space 
onto planes of atomic internal variables and momentum, we have found typical structures in 
Hamiltonian chaotic systems - - chains of resonant islands of different sizes imbedded in a 
stochastic sea, stochastic layers, bifurcations, and so on. The phenomenon of sticking of atomic 
trajectories to boundaries of regular islands found in numerical experiments should have a great 
influence to atomic transport in optical lattices. 

One of the aims of this paper was to describe analytically and numerically fundamental 
aspects of nonlinear dynamics of the atom-field interaction. We have done that to some ex- 
tent at the cost of simplifying the model. To be more realistic we should take into account 
spontaneous emission events. In this case, however, the equations of motion cease to be a deter- 
ministic dynamical system because they would include random terms. Our previous results on 
Monte Carlo modelling Hamilton-Schrodinger equations have shown much more complicated 
type of atomic motion which, besides of chaotic motion, caused by the fundamental atom-field 
interaction, includes a purely stochastic component caused by random events of spontaneous 
emission. We plan to study the effects of spontaneous emission on chaotic atomic motion in 
the future. 
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